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Correlation functions for the 2D random bonds Potts Models. 
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We study the spin-spin and energy-energy correlation functions for the 2D Ising and 3-states Potts model 
with random bonds at the critical point. The procedure employed is the renormalisation group approach of the 
perturbation series around the conformal field theories representing the pure models. For the Ising model, we 
obtain a crossover in the amplitude for the correlation functions which doesn't change the critical exponent. For 
the 3-state Potts model, we found a shift in the critical exponent produced by randomness. A comparison with 
numerical data is discussed briefly. 



1. INTRODUCTION 

For models with random bonds, the main prob- 
lem is to determine if the randomness leaves un- 
changed the critical properties of the pure sys- 
tem or if the singularities of the thermodynam- 
ical functions are modified. First results, which 
suggest an intermediate situation, have been ob- 
tained by Harris and Lubensk y ll[| , Grinstein and 
Luther ^ and Khmelnitskii |i| using the stan- 
dard 4 theory. Other cases, like long-range cor- 
related quenched defects have also been consid- 
ered, (see for example [|]||.) A first step in the 
understanding of the relevance of randomness was 
given by the Harris criterion |J : the randomness 
is relevant (irrelevant) if the specific heat expo- 
nent of the pure model is positive (negative) . Two 
dimensional systems are particularly interesting 
because of the rich structure of conformal invari- 
ance in this dimension. Assuming that a random 
model has a critical point with second order phase 
transition, the main interesting problem is to de- 
termine which conformal field theory represents 
this model at the infrared fixed point. Let us 
also mention that an exact result has been ob- 
tained by McCoy and Wu who considered a 
two-dimensional Ising model where only vertical 
bonds on a square lattice were allowed to acquire 
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the same random value. They found that the 
logarithmic singularity of the specific heat disap- 
peared completely. 

The models that we will study in this paper are 
the two dimensional Ising and Potts models with 
random bonds. For the case of Ising model, Harris 
criterion doesn't provide a qualitative answer of 
the relevance of the randomness (the specific heat 
exponent for the 2 — D Ising model is 0). First 
results for this model were obtained by Dotsenko 
and Dotsenko Q who showed that near the crit- 
ical point, this model can be represented by an 
n = Gross-Neveu model Q. With this tech- 
nique, they found that the specific heat singular- 
ity get smoothed as ln(ln(-^)) where \t\ is the re- 
duced temperature. Calculation of spin-spin cor- 
relation function by this technique which involves 
non-local fermionic representation of < a a >, 
was later questioned by Shalaev and Shankar 
fl^ , |ll| who gave arguments that the asymptotic 
behavior of this correlation function is unchanged 
by the randomness (see also Jl2[.) Let's also men- 
tion some other approaches. In pUfl , the situa- 
tion was questioned by Ziegler who claims that 
non-perturbative effects introduce an intermedi- 
ate phase around the critical point of the pure 
model. In an other approach, Mussardo and Si- 
monetti studied the exact formulation of the ran- 
dom bond Ising model in term of S'-matrix juj. 
On the other hand, numerical simulations of the 
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Ising model |l6|,[15| seem to confirm the theoreti- 
cal predictions of the specific heat and spin-spin 
correlation function asymptotic behavior. 

For the Potts model, the situation is more clear. 
Hear, the Harris criterion makes some precise pre- 
dictions: randomness is relevant and changes the 
critical behavior. Using conformal field theory 
techniques, Ludwig Jl7|] also perturbatively com- 
puted a shift in the critical exponent of the energy 
operator in the case of the random Potts model 
(this critical exponent for the Ising model is also 
unchanged) . 

The paper is organized as follows. In section 2 
we introduce the model. In section 3, we briefly 
present the methods for computing correlation 
functions, referring to Q for more details. Re- 
sults are presented in section 4 while in section 5 
a discussion of the results is presented and a com- 
parison with numerical data obtained recently. 

2. THE MODEL 

The partition function for a q-state Potts model 
with a fixed configuration of disorder is given 
by 



z(Jij) = E e 



(1) 



Here, Oi corresponds to the value of the spin at 
the location i on the lattice and takes the values 



o-j = 0, • • • , q - 1 



(2) 



and Jij corresponds to the coupling constants be- 
tween neighboring spins. These coupling constant 
can be separated in two parts 



Jij — Jo + 



(3) 



Jo is the coupling constant of the pure model 
(without disorder) while is the random part 
of the couplings. This partition function can also 
be written under the form 

Z(Jij) = 2 e - SoHS X)<*,i> s-i-i^ ( 4 ) 
{"<} 

where we have explicitly separate the non ran- 
dom part (to which corresponds a action asso- 
ciated to a conformal field theory at the critical 



temperature) and a random part. Thus, using the 
fact that the 5 rTi . r7j corresponds, in the continuum 
limit of the pure Potts model, to the energy field 
e(z), we can also express the partition function 
like 



Z(m) — 1 r{ a .ye J v ' ' 



(5) 



and m(z) is the continuum limit of A,j. Under 
this form we recognize a conformal field theory 
with a perturbation term J d 2 z m(z)e(z). 

The next step consists in performing the aver- 
age over the disorder i.e. over m(z). In fact, we 
need to compute the average over the free energy, 
i.e. over In Z(m). This is done easily by using the 
replica method. Taking the partition function of 
n identical copies of the system and analytically 
continuing to the limit n — > gives the quenched 
free energy 



-j3F = ln(Z) = lim 



Z n - 1 



n-*0 n 



where: 



n - } j S , a - / m(z) 2 t s a (z)d z 

Z n = ]]Tr {a ^ } e »=i J »=i ( 6 ) 

a=l 

the average of Z n is made with a Gaussian dis- 
tribution for m(z): 



Z^= I ' Y[dm(z)Z n e-^k {m(z) - mo)2 

J z 

which gives: 

n 

Z^ = l[Tr a . Sz e-S 

a=l 

and 

n p n 

S = 5Z^°' a + - 9 ° / £ a{z)sb(z)d 2 z 



(7) 



a=l 



m 



a=l 



a, 6=1 

e a (z)d 2 z (8) 



The terms in J 53 £ a(z)sb{z)d 2 z containing the 



.,6=1 



same replica label produce irrelevant operators 
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and can be omitted. Then S is reduced to 



S 



a=l 



TO 



So, a + 30 



e a (z)eb(z)d 2 z 



/n 



(9) 



In the limit too — > 0, this model corresponds 
to a conformal field theory perturbed by a term 
quadratic in the e operator. Then the "evolu- 
tion" of the coupling constants go and mo under a 
renormalisation group (R.G.) transformation can 
be analyzed as well as the behavior of the corre- 
lation functions. In the calculation of correlation 
functions < O(0)O(R) >, where O is some local 
operator, we will proceed perturbatively: 



< 0(0)0(7?) > 
< S7O(0)O(i?) > 



< O(0)O(R) > (10) 
~ < SfO(0)O(R) > +... 



where <>o means the expectation value taken 
with respect to Sq and 

Si = J H I (z)d 2 z = g Jj2 £ a(z)e b (z)d 2 z (11) 

The operator O is then renormalised as 

O -» 0(1 + A l9o + A 2 g 2 + A 3 g 3 + ■■■) = Z O 

The integrals of correlation functions involved in 
the calculation can be performed by analytic con- 
tinuation with the Coulomb-gas representation of 
a conformal field theory where the central 
charge is c = h + e'. The d term corresponds 
to a short distance regulator for the integrals. In 
addition, we also used an infrared (I.R.) cut-off r. 
The result is then expressed as an e' series with 
coefficients depending on r. The limit e' — > 
corresponds to the pure Ising model at the crit- 
ical point while the Potts model is obtained for 
some finite value of e'. We recall here some no- 
tations of the Coulomb-gas representation for the 
vertex operators 19 1. The central charge c will 
be characterized in the following by the parame- 
$ + e with 



2p 



c= l-24ao ; a± = a Q ± \/ 0% + 1 (12) 

= —1 



Note that for the pure 2D Ising model a\ = | and 
c = | while for the 3-state Potts model a+ = |, 
c — $ and e = —%r. 
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defined by 



The vertex operators are 



^nm(x) ^ 



(13) 



where 4>{x) is a free scalar field and where the 
a nm are given by 



n)a_ + — (1 — m)a + 



(14) 



The conformal dimension of an operator V nm (x) 



is /S. nm — Q-nm^nm with 

= 2a -a 



'-ran 



- (l+n)a_ + i (l+ro)a+(15) 



The spin field a can be represented by the vertex 
operator V^ iP _i whereas Vi,2 corresponds to the 
energy operator e. In the same way, we associate 
e ia + <t>(x) £ k e screening charge operator V+(x). 
Note that in the Ising case the a operator can also 
be represented by the V21 operator (since both 
operators coincide in the limit e — > 0). So, we 
can represent our spin operator by Vk,k—i where 
k = We have A = 2 for V 21 and A = ± for 

Vt).d-1 • 



3. Renormalisation Group Equations 

In this section, we will deal with the compu- 
tation of correlation functions of operators e and 
a. To compute them, one needs to determine the 
effect of the random coupling on the operators e 
and a and compute the renormalised operators e' 
and a' . This means that we want to compute the 
functions Z e and Z a such that 



e = Z e e and a' = Z^a 



(16) 



A convenient way to define Z e and Z„ is to con- 
sider the more general action 



n „ n 

a=l J o,6=l 



z)d 2 z 



a=l 



+ ///,, / > ]e a (z)d 2 z - h I ^2<r a (z)d 2 z (17) 

a=l 



This merely corresponds to the action used 
in (0) with an additional coupling of the a 
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field. Then, with the help of the operator al- 
gebra (O.A.) coming from contractions between 
e and a operators, we will compute the ef- 

n 

feet of the go J 53 e a (z)eb(z)d 2 z term (that we 

a, 6=1 

will denote go{ee) in the following for simplic- 

n 

ity) on the coupling terms mo / ^^£ a {z)d 2 z and 

0=1 

n 

/lo J ^3 a a (z)d 2 z. More precisely, we will com- 



pute 



a=l 



53 (go(ee)) l m Q / ^e a (z)d 2 z 
» \ J 0=1 y 

« n 



(18) 



and 



^(We))^o | 5>(*¥ 2 ^ 

-*y 53 ^w^ 2 



(19) 



m and ft. being the renormalised coupling con- 
stants. Obviously, this computation will be per- 
turbatively made only up to some finite power in 
<7o- In fact, the first step of the computation will 
be to determine the renormalised g constant on 
which Z e and Z a depend. 

3.1. Renormalisation of the coupling con- 
stant g. 

The renormalisation of the coupling constant 
g will be determined directly by a perturbative 
computation, g is also given by the O.A. produc- 
ing 

1 2 1 3 

9o (ee) + 2 (3o (ee)) + g (j9o(ee)) +'•• 



P n 

1 / ^2 e a{z)e b {z)d 2 z 



,,6=1 

_2 1 A „3 



with <7 = go + ^2<?o + ^3fo + ' ■ ' where A2 comes 
from 

2 / 53 M£)£&(2)d 2 2 / 53 £c(^)£d(^)d 2 z 

** 1 ** 1 



/n 
53 £a(2)ei(«)d a «+---(20) 
„ 1 



a, 6=1 



and A3 from 



l ( / 53 £ a(2)£fc(z)d 2 Z I 
\ a, 6=1 / 

/n 
^ £ Q (z)£ b (z)d 2 z + --- (21) 

a,6=l 

The computation of A2 is made by contracting 
two e operators. We obtain (with b = c ^ a,d) 

~ n p n 

5 / 53 £ a(a;)£6(a;)d 2 2: / 53 £ c(y) £ d{y)d 2 y 

*' I, 1 J „ j 1 



2(n-2) / <£(^)£(y) > ^x 

|x— 3/ 1 <r 



/ 53 e Q (a;)e d (a;)d 2 a; (22) 



4?r(n - 2) 



,,d=l 



/ ^rfc / 53 £a(^£d(^)rf ; 

J y< rV J a,d=l 

-3e 



-47r(n - 2)- 



3e 



53 £ a (a;)£<i(a;)d 2 



a.d=l 



Thus A 2 = — Am(n — 2)^^^— . Computations at 
higher order will become very complicated. In 
the following, we will just present the results, re- 
ferring to [jl8j for all the technical details. For the 
computation of g(r), calculations were made up 
to the third order with the following result : 



g(r) = r -^(g Q -g 2 47r(n-2) — 

,„ o. ^r~ 6e / 2(n-2) 
+ gl%ir 2 {n-2)— I 1 



(23) 



3e 



3e 



Note that we multiply the result by r~ 3e in order 
to obtain a dimensionless coupling constant g(r). 

From there, we can compute directly the /3- 
function: 



Kg) 



dg 



-3eg(r) + Air(n - 2)g 2 (r) 



a, 6=1 



c,d=l 



dln(r) 

- 16ir 2 (n-2)g 3 {r) + 0(g i (r)) (24) 
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Finally, taking the limit n — > 0, we obtain for the 
/3-function up to the third order : 

p(g) = -3eg - Sng 2 + 32tt 2 9 3 (25) 

From this /^-function, we determine immediately 
the effect of the disorder on the model. For the 
Ising model (e = 0), the infrared fixed point is 
g = 0, while for the 3-state Potts model (e < 0), 
a new infrared fixed point is reached with g c = 

_ 11 + 9gi + 0(e 3) 

3.2. Renormalisation of a and e 

In order to be able to compute the correlation 
functions of a and e, the second step is to de- 
termine the effect of the renormalisation on these 
operators. One needs to compute the multiplica- 
tive functions Z a and Z z . This will be made 
by computing the renormalised coupling terms 

n p n 

to J ^^e a {z)d 2 z ~ (moZ £ ) / ^^e a {z)d 2 z and 

a=l ^ a=l 

n r, n 

h IzZ a ^ d2z = / zZ a ^ d2z - A di " 

a=l J a=l 

rect computation of both to and h will provide us 
with the functions Z a and Z e . As for the compu- 
tation of g = Z g go, we will compute in perturba- 
tion : 



/n „ n 

^2e a (z)d 2 z + (g (ee)m / ^e a (z)rf 2 z 
a=l J a—1 

1 f ™ 

+ ^(3o(££)) 2 "7 / ^ e a {z)d 2 z + ■■■ 

1 J 0=1 

/n 
^e Q (z)d 2 z (26) 

a=l 

and to = m (l + -B1.90 + B 2 g 2 + • • •) with B\ 
defined by 

/n - n 

^2 £ a (z)s b {z)d 2 z / ^e a {z)d 2 z 
- •— 1 ^ a=l 

„ n 

= B x j J> a (z)d 2 z (27) 



a, 6=1 



and B2 by 



/ £ a (z)e h (z) / ^2e a (z)d 2 z 

1 a, 6=1 / a—1 



= S 2 / ^e Q (^ 2 z (28) 
0=1 

The details of the computation are presented in 
fllgl, with the result 



r +2 to 



(r) = 7770 ^1 — 47t(t7 — 1) 



-3e 



.90- 



3c 



r" 6e , 477 - 6. 



+4^(77 -l)g 2 — (l + ——) (29) 



3e 



3c 



Here again, we multiply m(r) by r~ 1+ ^ e in order 
to obtain a dimensionless coupling constant. The 
R.G. equation for Z E is thus given by 



dln(Z £ (r)) 
dln(r) 



4ir(n-l)g-8ir 2 {n-l)g 2 (30) 



Similarly, for the coupling constant h$, we com- 
pute up to the third order : 

^ p n \ „ n 

/ ^2 £ a( z )£b{z) / ^(J a (z)d 2 Z 
. a, 6=1 / a — 1 

« n 

= Q j J2°a(z)d 2 Z (31) 



We give here directly the result : 



-6c 



r"r-»W&( r ) = ^(1 + (77 - l) ffo V — x (32) 

4 r 2 f-2)r 2 fi) 



-12(77 -1) (77 -2)c^ 3 



-9e 



9c 



( 1 + 9( 2 - A ) r 2(-i)r2(-i) )) 

The multiplicative term r~^~ a ( e ) in front of h(r) 
is introduced in order to make this parameter di- 
mensionless. Here, a(e) is a function of e depend- 
ing on which representation of the spin field we 
are taking in the Coulomb gas picture (see sec- 
tion 2). Its explicit form will be irrelevant in the 
following. The corresponding R.G. equation for 
Z a will be given by 



dln{Z a (r)) 
dln(r) 



(33) 



G 



-3(n- l)g 2 (r)n 2 e 



i + 4 (2 \) r2 H)r 2 (§) 

+4(n- l)(n-2)7rVW 
4. Correlation Functions 

We now have all the ingredients needed in order 
to compute the correlation functions. They will 
be calculated with the help of the R.G. equations, 
for the theory with mo, ho — > 0. From the R.G. 
equations, we have : 



< s(0)e(sR) > r ,g(r) = 

Z 2 (sr,g(sr)) 



(34) 



-2A e 



Z 2 {r,g(r)) 
This can be written as 



< £(0)e(R) > r ,g(sr) 



< e(0)e(si?) > r , 9 (r) = 



(35) 



2 / ^ 



e «> s 2A - < e(0)s(R) > r , g ( sr ) 

where we used the notation : 



dlnZ £ 
dlnr 



(36) 



and g(s) = g(sr);go — g(r). We assume now r 
to be a lattice cut-off scale. In a similar way for 
< a(0)a(R) > the R.G. equation is : 



< a(0)a(sR) > r ,g( r y 



(37) 



2 / 



-2A a 



< cr(0)cr(.R) > r , s ( sr ) 



with 
dlnZfj 
dlnr 



7<r(ff) 



(38) 



In equations (|3q)-(|37|), i? is an arbitrary scale 
which can be fixed to one lattice spacing r of 
a true statistical model. The dependence of 
< cr(0)a(r) > r ,g(s) on s wn l then be negligi- 
ble, assuming that there are not interactions on 
distances smaller than r. Therefore, it reduces 
to a constant. Then, s will measure the num- 
ber of lattice spacings between two spins in < 
a(0)a(sR) >. In the following, we adopt the 
choice r = 1. 



4.1. The Ising model 

The Ising model corresponds to the case e — ► 
and so the /3 function is : 



Pig) 



2„3 



-8irg z + 32tt z s 



(39) 



Therefore, we can see that the I.R. fixed point is 



located at g = 0. Also we have, by eqs_. 
for n = 0, e = and definitions (J3q) , 



7e(fl0 

7a (5) 



2„2 



-47r<7 + 87r g 



B5D. © 

(40) 
(41) 



The integral for the e correlation function, 



1 — 2irg dg 



eq.(f35[), gives : 




ffW 


9« 




■ /" 




go 








' /< 




90 



1 - 4tt^ 

/ (1 + 2^5) — 
30 

Mff(s) - g Q ) + ln 



(42) 



5(f) 
5o 



Now, we need to compute <?(s). The integration 
of equation /5(g) = gives : 



90) 



dg 



-8irg 2 + 32ttV 



dli 



with the following solution up to the second or- 
der : 



g{s) 



go 



(43) 

^ + 47rffo/n(l + 87Tffo/n(s)) ^ + Q ^ 



1 + 8Trgoln(s) I + 8irgoln(s) 

So, < se > correlation function is given by : 

g(s) 



< e(Q)e(«) >9o 
I 



50 



(1 - 27r( 50 - 5(s))) s~ 2A 



x (44) 



1 + 8ngoln(s) 1 + 8ngoln(s) 

(ln(l + &TTg ln(s)) - 4Trg ln(s))) + 0(g 2 ) 

For the a correlation function, the computation 
is similar. In fact, in that case, we have J a (g) = 



7 



8ir 3 g 3 . Thus, keeping only the first order of the 
/3-function (i.e. (3(g) = —8irg 2 ), we obtain : 



90) 



dg 



90) 
-2tt 2 J gdg 

9o 



= -* 2 (g{s) 2 -gl) + 0(gl) 



(45) 



The < (7cr > correlation function is then found to 
be given by : 



<a(0)a(s) >~ (1 + tt 2 (g 2 - g(s) 2 )) s 

1 



1 + 7r 2 .9o 



1 - 



(1 + 8irg a ln(s))< 



-2A„ 



,-2A„ 



(46) 
O(5o 3 ) 



The calculation of the g 3 term in the (3 function 
and the g 2 term in the renormalisation of e was 
already done in jl^] , extending the one loop result 
of Q . We recovered these higher order corrections 
using a different technique, which allowed us to 
calculate also the modified correlation function of 
the spin operators. 

4.2. The Potts model 

We consider here the 3-state Potts model. 
With our conventions this case corresponds to 
e = — j^. (3(g) is given in eq.(|||), and 



ie(g) 
lAg) 



-Airg + 8n 2 g 2 



l + 3 (2 X \ 2 (-\)Y 2 (-\), 

3.3 (48) 



(47) 
9 2 



8^g 



At long distances, the integrals in eqs.(|35|), ( |37[ ) 
will be dominated by the region g ~ g c , with 
g c = + ff^ + 0(e 3 ). This is different from 
the Ising model, because here j e (g c ) and 7c- (g c ) 
have finite values for g — g c . Thus, 



sO) 



ie(g) 
(3(g) 



dg ~ j £ (g c )ln(s) 



(49) 



90 

and 

90) 

/ 

9o 



i<y(g) 
(3(g) 



dg ~ le(g c )ln(s) 



(50) 



The correlation functions can then be deduced 
directly : 



< e(0)e(a) >~ s -( 2 ^- 2 ^)) 
and 

< a(Q)a(s) >~ s -(2A.-2 7 .( ffo )) 



(51) 



(52) 



So, we can see that a direct consequence of the 
new IR fixed point is a modification of the critical 
exponents A e and . A straightforward compu- 
tation will give these new exponents : 

2A^ = 2A £ -2 7e ( 5c ) 

= 2A £ + 8ng c - 16tt 2 5 2 



2A £ -3e+^e 2 + (9(e 3 ) 



(53) 



and 



2A ff = 2 Act - 2 la (g c ) 
= 2A a - 6ir 2 g 2 e(l 



(54) 



r 2 (-f)r 2 (i) 



Here A = \ and the final result for the new critical 
exponent is : 



2A' = 2A„ 



27 r 2 (-|)r 2 (i) 
i 6 r 2 (-i)r2(-i)' 



0(e 4 ) (55) 



The value of the critical exponent for the energy 
operator was already computed, up to the second 
order, by Ludwig [jl7|]. While for the pure case 
A e = 0.8, here we obtained A.' s = 1.02 + 0(e 3 ). 
The difference between these exponents is quite 
important and should be measured in a numeri- 
cal simulation of this model. See next section for 
a discussion on this point. For the spin opera- 
tor, deviation of the critical exponent from the 
pure case appears only at the third order, which 
we compute. This deviation is in fact very small. 
While in case of 3-state Potts model without dis- 
order 2 A CT = y|r, we obtain the new critical expo- 
nent 2 A^ = jl+O; 00264 = 0,26931. Thus, the 
deviation corresponds to a modification of 1%. 
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5. Discussion 

In case of Ising model spin-spin function the 
calculation of up to third order of the renor- 
malisation group was needed to find the devia- 
tion from the perfect model case, in the form of 
the cross-over in the amplitude. This completes 
the observations of (l(]-[l2|, based on absence of 
renormalisation of this function in the first order, 
that asymptotically the spin-spin function has the 
same exponent as in case of the perfect lattice 
model. 

Recently the numerical simulations of the ran- 
dom Ising model has been performed which mea- 
sure directly the deviation of < aa > from the 
pure Ising model at the critical point |fuj| . These 
measurements were made for disorder such that 
87rgo ~ 0.3. Deviations predicted by our compu- 
tations are very small. They correspond to 0.1%. 
The deviations obtained in numerical simulations 
are around ten times larger, and they are of oppo- 
site sign, i.e. correspond to an extra decrease of 



the spin-spin function with distance r. In |16|, it 
has been checked that this decrease corresponds, 
within the accuracy of the measurements, to a 
factor function of the ratio r/L, F(r/L), r being 
the distance between the spins and L is the lattice 
size. So they correspond to finite size effects, be- 
ing different for perfect and random models. We 
would suggest, on the bases of our calculation of 
the r dependence of the spin-spin function on an 
infinite lattice, that numerical deviations will con- 
tinue to be plotted by the same curve F(r/L), if 
one measures < era > for different lattice sizes as 
it has been done in [jl6| , until the accuracy reaches 
the value of the r-deviation which we calculated 
here. Only then the curves for different L will 
split. 

Recently numerical simulations of the 3-state 
Potts model with disorder were performed by one 
of the authors pfj. In these simulations, A a > was 
measured as well as A e /. For A a > , it is very 
difficult to obtain any conclusion. The measured 
value A a > is very close to A a and error bars do 
include the case A a > = A a up to a deviation of 
1%, which is the deviation predicted by analyti- 
cal computations presented in this work. For the 
energy-energy correlation function, the measured 



exponent is A £ r = 1.065 ± 0.02, to be compared 
with the prediction of our computations: A £ i = 
1.02. This is quite a satisfactory result if we re- 
member that the result A E i = 1.02 is a perturba- 
tive computation while the pure model (around 
which we perturbed) has the value A e = 0.8 
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